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A NOTE ON THE FRACTIONAL LOGISTIC EQUATION 


I. AREA, J. LOSADA, AND J. J. NIETO 


Abstract. In this short note, we show that the real function recently proposed by 
Bruce J. West [Exact solution to fractional logistic equation, Physica A 429 (2015) 
103-108] is not an exact solution for the fractional logistic equation. 


1. Introduction 

The topic of Fractional Calculus (that is, the study of integrals and derivatives of non¬ 
integer order) is nowadays a growing interest area of mathematics. Moreover, physicists 
and engineers are also really interested in applications coming from this nice theory 
which origin goes back to the theory of differential calculus; more precisely to Leibniz’s 
famous note in his letter to L’Hopital dated 30 September 1695. 

For a long time, the theory of Fractional Calculus developed only as a theoretical 
held of mathematics. However, in the last decades, it was shown that some fractional 
operators describe in a better way some complex physical phenomena, especially when 
dealing with memory processes or viscoelastic and viscoplastic materials. Well known 
references about the application of fractional operators in rheology modeling are nna. 
One of the most important advantage of fractional order models in comparison with 
integer order ones is that fractional integrals and derivatives are a powerful tool for the 
description of memory and hereditary properties of some materials. Notice that integer 
order derivatives are local operators, but the fractional order derivative of a function in a 
point depends on the past values of such function. This features motivated the successful 
use of fractional calculus in CRONE im and PID controllers [T^ Chapter 9]. 

Due to this, in the last two or three decades, a great interest has been devoted to the 
study of fractional differential equations. Thus, fractional order differential equations 
(one day called extraordinary differential equations) play today a very important role 
describing some real world phenomena. For a complete exposition of the theory of 
Fractional Calculus one can see [3 El ini [n]. 

The exponential function, exp(f), plays a fundamental role in mathematics and it is 
really useful in the theory of integer order differential equations. In the case of fractional 
order, it loses some beautiful properties and Mittag-Leffler function appears as its natural 
substitute. Next, we recall its dehnition and some basic properties. 

Definition 1. The function Ea{z) is named after the great Swedish mathematician Gosta 
Mittag-Leffler (18f 6-1927) who defined it as a power series given by 

OO ^ 

= E fuuu)' “ > 
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This function provides a simple generalization of the exponential function because of 
the replacement of k\ =r{k + 1) by {ak)\ = r{ak + 1) in the denominator of the terms 
of the exponential series. Due to this, such function can be considered the simplest 
nontrivial generalization of exponential function. 

Definition 2. A two parameter function of Mittag-Leffter type is defined by 

OO ^ 

= E T{ak + py “ > 0. -3 > O' 

During the hrst half of the twntieth century, the Mittag-Leffler functions remained 
almost unkonown to the majority of the mathematicians. However, recently attention 
of mathematicians and other scientists towards functions of Mittag-Leffler type has in¬ 
creased and today some mathematicians like to refer the classical Mittag-Leffler function 
as the Queen Function of Fractional Calculus, and to consider all related functions as 
her court. 

Recently, a complete monograph (see 0 ) about Mittag-Leffler functions and its several 
applications has been published. More information about this kind of special functions 
ans its relations with Fractional Calculus can be found in [S]. 

It follows from previous definitions that, for instance. 


Ei,i{ 


= exp ( 2 ;), 

E2,i{z^) 

= cosh(^). 
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) = cos(z). 

» 

9 ’ 

E2,2{z"‘ 


which proves the importance of Mittag-Leffler functions in mathematics. Furthermore, 
using the following well known result 

we can easily see that with A G M, is a nontrivial eigenfunction for Caputo 

fractional derivative operator (previously denoted by see [71 Section 2.4]). That is, 
if f{t) = F^a(Af") for t > 0 then 

°D“/(t) = A/(t), t>0. 

In the past, some authors used to derived some results the following property of Mittag- 
Leffler function 

^a(a(t + s)“) = E„(an^a(as“), 

where a is a real constant. Unfortunately, such property is unavailable unless a = 1 or 
a = 0, both of them are trivial situations, since we obtain exponential function in the 
first case, and a constant function secondly. For more information about the invalidity 
of such equality we refer to HD. 

In this note we discuss a similar property related again with Mittag-Leffler function 
and using such property, we conclude that the function introduced in [TB] is not an exact 
solution for the fractional logistic equation. It may be obvious, but nevertheless Mittag- 
Leffler function shares some analogue properties with exponential function, these kind of 
results show that exponential function cannot be replaced, in general, by Mittag-Leffler 
function in the theory of Fractional Calculus. 
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2. Main results 

In 1838, P. F. Verhulst introduced a nonlinear term into the rate equation; he was 
studying population models and he wanted to avoid the catastrophic predictions previ¬ 
ously proposed by T. Malthus, who had used the rate equation to model human pop¬ 
ulation growth (for historical references see [IS]). By this way, P. F. Verhulst obtained 
what today is known as the logistic equation: 

(1) u'{t) = ku{t){l — u{t)), t>0. 

This differential equation is one of the few nonlinear differential equations that has a 
known exact closed form solution, which is given by 

'w(f) = -- ~\ - / —TT) ^ > 0, 

Mo + (1 - uo) exp(-fct) 

where uq is the initial state (that is, m(0) = uq = iV(0)/iVmax, where iV(0) is the total 
population at the initial time and iVmax is the carrying capacity of the ecosystem). Since 
that time, the logistic equation has found several applications ranging from machine 
learning to modeling growth of tumors. 

Fractional logistic equation has been studied, see for instance but until now no 

exact is known. Recently, Bruce J. West has published a research paper (see US]) where, 
using the Carleman embeding technique, he proposes an exact solution for the fractional 
logistic equation. More concretely, in such reference it is said that the function given by 


(3) 


u{t) = 


n=0 


Uq-I 

Uq 


Ea,{-nk^r), t > 0, 


may be an exact solution of the fractional differential equation 

(4) (t) = k°'u{t)(l — u{t)), t >0, 0 < a < 1. 

Of course, following [16], notice that if a = 1 then function given in ([^ is equivalent 


to 


u{t) = 


n=0 


Up - I 
Up 


exp(—n/cf) = 


Up 


t > 0 


Mo + (1 - "lio) exp(-H) ’ 

which coincides with (|^. That is, if in (|^ a = 1, then we obtain the solution of the 
integer order logistic equation ([^. Moreover, fractional solutions given by (|^ have 
similar properties to integer order solution given by ([^. 


Remark 1. Since Carleman embedding technique was introduced to solve integer order 
differential equations, not fractional order differential equations, the author of [TB] warms 
us about the necessity to validate solution of 0 given by 0. This is the main objective 
of this short note. 


In what follows, we assume that function given by (|^ is an exact solution of the 
fractional differential equation (|^. 

Lemma 1. If the function given by ^ is an exact solution of fractional logistic equation 
0 of order 0 < a < 1, then for all n E N U {0} we have that 

n 

(n + 1)E„(-mA:“s“) = E„(-(n - j>“f“)E„(-jrt“). 

j=p 


( 5 ) 
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Proof. Consider for each j G N U {0} 


Qj — 


Up - 1 
Uq 


Eo^i-jk^n. 


Thus, function given by ([^ can be written as 

OO 

u{t) = On, t > 0. 


71=0 


Moreover, given 0 < a < 1, we have that 






n=0 


Uo 


(t) = -k^^[ 


n=0 


Uq 


and using the Cauchy product, we obtain: 
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it)il-uit)) = J2('^^) E^i-nk'^r) 


n=0 

OO 


Uo 


n=0 

OO 

z 

n=0 


Uo 

Up-I 
Uo 


71=0 
OO n 


5 :(^) E^(-nk‘e)-Y.j: 


Uo 




71=0 j=0 

71 

E^{-nk^n -j2Ea{-{n- j)k^nE^{-jk^n 


3=0 


Since u{t) is a solution of equation (|^, we deduce that for all n e N U {0} it must be 

71 

(n + ^ ^a(-(n - 

3=0 

which is the equality given in the statement of the lemma. □ □ 

Proposition 1. In the condition of Lemma\^ we have that 
(6) E^ (-2fc“t“) = E^ i-k^r) i-k^E). 


Proof. Considering n = 2 in equality ([^, we obtain the desired result. 
Proposition 2. Identity ^ of Proposition^ holds if and only if a = 1. 

Proof. By definition of the Mittag-Leffler function we have that, 

-2k^t^T 


□ □ 


Ej-2k^E) = 

71=0 


^ ^ /-r(na + l)’ 


(7) 


t > 0. 
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Moreover, using the Cauchy product, we obtain: 




( 8 ) 


+ i)r(ja + 1) 




n=0 

oo 


j=0 

n 




r((n - j)a + i)r(jQ; +1) 

2 -n 


n=0 

oo 


j '^0 r((n-j> + l)r(ja + 1) 


t > 0, 


n=0 


where 

^ ^ 2~”' 

§r((n-j> + l)r(ja + l)‘ 

Equality ([^ implies that coefficents of series expansions Q and ([^ must be equal. 
Particularly, if we choose n = 2, we must have 


1/2 1 \ _ 1 

4 Vr(2a + 1) ^ r(a + l)r(a + 1 ) y r(2a + 1) ’ 


or equivalently. 


(9) r(2Q; + 1) ^ 1 

^ ^ 4r(a + l)r(a + l) 2‘ 

Since if a G (0,1) then r(2a + 1) < r(3) = 2 and 4r(a + 1)^ > 4r(l)^ = 4, it is clear 
that equality in ([^ holds if and only if a = 1. □ □ 

Remark 2. From Proposition it is dear that function given by is a solution of 
differential equation 0 if and only if a = 1, but in such case, we have an integer 
order differential equation and function defined by coincides with the solution of such 
differential equation. 

Notice also that since when a = 1 then Ea(t) = exp(f), identities @1 and ^ are well 
known properties. 


Remark 3. Differentiating in both sides of equation 0, we deduce another identity 
which only holds if a = 1, 

Ea,a (-2fc“r) = E^ {-k^r). 


Using Mathematica HZ] (see also m), we can reaffirm analytical results previously 
obtained. In hgure 1 some graphs can be found in this direction. Notice that when a is 
near 1 the difference is smaller than in the case a near 0. 


Recently, Turalska and West US] have used a spectral technique to study the solution 
to the fractional logistic differential equation since the exact solution is not yet known. 
For the approximate siolution there is a small, apparently generic, systematic error that 
they are not able to fully interpret. 
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alpha=0.9 alpha=0.75 



Figure 1. Comparison considering diferent values of a {a = 

0.9, 0.75, 0.5, 0.25) between Ea (—(continuous line) and the prod¬ 
uct Ea Ea (dashed line). 

3. Conclusion 

In this short note, we have shown that the function recently proposed by Bruce J. 
West is not an exact solution for the fractional logistic equation. 

In any case, we would like to mention that the use of Carleman embeding technique 
may be a powerful tool in the study of some fractional differential equations, but its 
application must be studied carefully. 
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